Analytical solution of the equations describing interstitial migration of impurity atoms 
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An analytical solution of the equations describing impurity diffusion due to the migration of 
nonequilibrium impurity interstitial atoms was obtained for the case of the Robin boundary con- 
dition on the surface of a semiconductor. The solution obtained can be useful for verification of 
approximate numerical solutions, for simulation of a number of processes of interstitial diffusion, 
and for modeling impurity diffusion in doped layers with the decananometer thickness because in 
these layers a disequilibrium between immobile substitutionally dissolved impurity atoms, migrating 
self-interstitials, and migrating interstitial impurity atoms can take place. To illustrate the latter 
cases, a model of nitrogen diffusion in gallium arsenide was developed and simulation of nitrogen 
redistribution from a doped epi-layer during thermal annealing of GaAs substrate was done. The 
calculated impurity concentration profile agrees well with experimental data. The fitting to the 
experimental profiles allowed us to derive the values of the parameters that describe interstitial 
impurity diffusion. 

PACS numbers: 07.05.Tp; 02.60.Cb; 61.72.Tt; 66.30.Dn; 66.30.Jt 

Keywords: diffusion modelling; equation solution; interstitial; nitrogen; gallium arsenide 



I. INTRODUCTION 



At present numerical methods have been widely used 
for simulation of solid state diffusion of ion-implanted 
dopants and impurity atoms added during epitaxial 
growth (see, for example, [3, [2, 3]). As a rule, to simulate 
the diffusion of impurity, a system of equations describing 
a coupled diffusion of different mobile species and their 
quasichemical reactions during annealing is solved. Due 
to a great number of differential equations and the com- 
plexity of the system as a whole, the problem of the cor- 
rectness of a numerical solution is very important. One 
of the best ways to verify the correctness of the approxi- 
mate numerical solutions is a comparison with the exact 
analytical solution of the boundary value problem under 
consideration. Such analytical solutions can be derived 
for the special cases of dopant or point defect diffusion. 
For example, in Ref. [1] an analytical solution for the 
point defect diffusion based on the Green function ap- 
proach was obtained. It was supposed in [i| that nonequi- 
librium point defects were continuously generated during 
ion implantation of impurity atoms and diffused to the 
surface and into the bulk of a semiconductor. The sur- 
face was considered to be a perfect sink for point de- 
fects. In Ref. [5] a process of impurity diffusion during 
ion implantation at elevated temperatures was investi- 
gated analytically. It was supposed that the implanta- 
tion temperature was too low to provide a traditional 
diffusion by the "dopant atom - point defect" pairs, but 
was enough for the diffusion of nonequilibrium intersti- 
tial impurity atoms to occur. Unlike [4], in Ref. [j| a 
system of equations, namely, the conservation law for 
substitutionally dissolved impurity atoms and the equa- 
tion of diffusion-recombination of nonequilibrium inter- 
stitial impurity atoms have been solved analytically by 
the Green function approach. Reflecting boundary con- 
dition at the surface of a semiconductor has been chosen 



to describe the interaction of interstitial impurity atoms 
with the interface. Due to this condition, a diffusion 
problem has become symmetric with respect to the point 
x = 0. For simplicity, the condition of zero impurity con- 
centration for x — ► ±oo has been used. It is interesting 
to note that analytical solutions for gold diffusion in sili- 
con due to the Frank- Turnbull and kick-out mechanisms 
were obtained in Ref. Q and Refs. @, 0] , respectively. It 
was supposed that there was a local equilibrium between 
substitutionally dissolved impurity atoms, vacancies (or 
sclf-interstitials for the kick-out mechanism) and intersti- 
tial impurity atoms. The case of nonequilibrium intersti- 
tial impurity atoms was not considered in these papers. 
The very interesting case of coupled diffusion of vacancies 
and self-interstitials was investigated in [1, Q . The equa- 
tions of the diffusion of vacancies or of sclf-interstitials 
are similar to the equation of diffusion of impurity in- 
terstitials. However, the solutions obtained in [H, [9[ are 
difficult to use for describing the impurity diffusion gov- 
erned by nonequilibrium interstitial impurity atoms, be- 
cause a condition of local equilibrium was used in all 
these papers. Besides, a generation rate was assumed to 
be equal to zero in [9( or equal to a constant value in 
Q. In Ref.[T3| analytical one-dimensional (ID) solutions 
of the equations that describe impurity diffusion due to 
migration of nonequilibrium impurity interstitials were 
obtained for the case of impurity redistribution during 
ion implantation at elevated temperatures and for diffu- 
sion from a doped epitaxial layer. The reflecting bound- 
ary condition on the surface of a semiconductor and the 
conditions of constant concentration on the surface were 
used in the first and second cases, respectively. These an- 
alytical solutions were obtained on the finite-length ID 
domain that is very convenient for comparison with nu- 
merical solution. Moreover, simulation of hydrogen diffu- 
sion in silicon during high-fluence low-energy deuterium 
implantation at a temperature of 250 °C and beryllium 
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diffusion from a doped epi-layer during rapid thermal an- 
nealing of InP /InGaAs heterostructures at a temperature 
of 900 °C was carried out on the basis of the analytical 
solutions obtained. It is interesting to note that the cal- 
culated impurity concentration profiles agree well with 
the experimental data that made it possible to derive 
the parameters of interstitial diffusion. This means that 
the solutions obtained are useful for solving a number of 
problems of solid state diffusion. The main goal of the 
present work is to continue the investigations of 0, [l(| 
to obtain a similar analytical solution for a more intri- 
cate case of Robin boundary condition on the surface of 
a semiconductor. 



II. ORIGINAL EQUATIONS 

As in Refs. [H, [To| , it is supposed that the processing 
temperature is too low to provide diffusion of substitu- 
tionally dissolved impurity atoms, but is high enough for 
the diffusion of impurity interstitials to occur. Nonequi- 
librium interstitial impurity atoms can appear due to ion 
implantation or due to the replacement of the impurity 
atom by self-interstitial from the substitutional position 
into the interstitial one (Watkins effect [HI), or due to 
dissolution of the clusters or extended defects which in- 
corporate impurity atoms. It is also supposed that the 
concentration of impurity in the doped regions originat- 
ing from the migration of nonequilibrium impurity in- 
terstitials is smaller or approximately equal to rij. Here 
rii is the intrinsic carrier concentration at the processing 
temperature. If the concentration of substitutionally dis- 
solved dopant atoms C is higher than rij, let us confine 
ourselves to the case of neutral impurity atoms in the 
interstitial position. Then, the system of equations de- 
scribing the evolution of impurity concentration profiles 
includes [El. Hoj: 

(i) the conservation law for substitutionally dissolved 
impurity atoms: 



8C(x,t) _ C AI (x,t) 



+ G AS (x,t), 



dt 
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(1) 



(ii) the equation of diffusion for nonequilibrium inter- 
stitial impurity atoms: 
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I A! = Vd^W^ , ~g AI (x,t) = G AI (x,t)r 
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Here C AI is the concentration of nonequilibrium im- 
purity interstitials; G AS is the rate of adding of impurity 
atoms, which immediately occupy the substitutional po- 
sitions, or (with the negative sign) the rate of the loss 
of substitutionally dissolved impurity atoms due to their 
transfer to the interstitial position; d AI and t ai are the 
diffusivity and average lifetime of nonequilibrium inter- 
stitial impurity atoms, respectively, and G AI is the gen- 
eration rate of interstitial impurity atoms. We use a 
steady-state diffusion equation for impurity interstitials, 
because of the relatively large average migration length of 
nonequilibrium interstitial impurity atoms (Iai 3> I fail, 
where lf a ii is the characteristic length of the decrease 
in the concentration of substitutionally dissolved dopant 
atoms) and due to the small average lifetime of nonequi- 
librium impurity interstitials tai in comparison with the 
duration of thermal treatment tp. 

The system of equations (JTJ) , © or (fl]) , © describes 
impurity diffusion due to migration of nonequilibrium in- 
terstitial impurity atoms. To solve this system of equa- 
tions, appropriate boundary conditions are needed. Let 
us consider, in contrast to 5], the finite-length one- 
dimensional (ID) domain [0, xb], i.e., the domain used 
in ID numerical modelling. The Robin boundary condi- 
tions on the surface of a semiconductor 
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is added to E q. (EH) in comparison with the solutions 
obtained in [f| [l(| • A similar condition in the bulk of a 
semiconductor 



R ,4/ C 

w"d A1 



AI 



dx 



x = x B 

as well as the initial conditions 



xb 



3 i 



(6) 



and w? , Wo , w? are the constant 



C(x, 0) = C (x) , C A1 {x,0) = C A1 = const (7) 

are also used for the final formulation of the boundary 
value problem. 

Here, wf, wf , w 3 anu w x , w 2 , w 3 
coefficients specifying the concrete type of real boundary 
conditions; C AI is the equilibrium value of concentration 
of interstitial impurity atoms in the bulk of a semicon- 
ductor (it is supposed that C A J is equal to zero for many 
cases of interstitial impurity diffusion) . 

To derive an analytical solution of this boundary value 
problem, the Green function approach [l2| can be used. 



III. ANALYTICAL METHOD AND SOLUTIONS 

The suggestion about the immobility of substitution- 
ally dissolved impurity atoms allows one to solve Eq. ((T|) 
independently of Eq. © or Eq. Q: 
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and on the right one: 



C{x,t) 



-AI 



C A1 {x,t)dt + G Ab {x,t)dt + C (x) 



(8) 

We will use expression ([5]) together with the steady- 
state solution of Eq. ^ obtained by the Green function 
approach 



xb 



C AI (x,t) = I G{x,S)w{t,t)dt, 



(9) 



Q2(xb) =w 1 d 



BjAI 
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(17) 



Following we can write the functions ws{x) and 
wb{x) as 



w s (x) 



<y(-i)wf if wf ^ . 



—gS (-z)wf ifwf^O. 



(18) 



where the standardizing function w(x,t) [12[ has the fol- 
lowing form: 



I 2 



+ w s (0 + w B (0- (10) 



Here G(x,£) is the Green function for Eq. ([3]). Using 
the standardizing function w(x, t) allows one to reduce 
the previous boundary value problem to the boundary 
value problem with zero boundary conditions: 



w B [x) 



wfd AI 



6(x B - x)wf if wf ^ . 



b<5 (x B - x)wf ifwf^O. 



(19) 



Let us consider the following Robin boundary condi- 
tion on the surface of a semiconductor [x = 0): 
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x = xb 



^C AI \x = x B =°- ( 12 ) 



The Green function for Eq. ([3]) with boundary condi- 
tions (|11|) and (fT2|) has the following form [T^ |: 
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where 
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K = Q 1 (x)Q 2 (x) - Qi(x)Q 2 (x) = const . (14) 
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i.e., wf = -1, wf = vf f f, wf 
' Dirichlet boundary condition 



0, and, for simplicity, 



C(x B ,t) = C A 



AI 



(21) 



in the bulk of the semiconductor, i.e., wf = 0, wf = 1, 
wf = C AI . Here, vf^ is the effective escape velocity 
specifying the intensity of impurity evaporation from the 
surface of the semiconductor. 

Then, the solutions Q\ and Q 2 have the form 



x 

Iai 



Qx(x) = -d AI cosh ( -Z- ) - l AI v s eff sinh [ f- ) , (22) 



x 
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Here Q\ and Q 2 are the linearly independent solutions 
of the homogeneous equation 
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with the following conditions on the left boundary: 



Qi(0) =w?d 
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= const , 
(24) 
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ws(x) = 0, w B {x) 



-S' (x B - x)Cf 



(26) 



Let us consider a buried layer highly doped with impu- 
rity atoms. Such a layer can be formed by ion implanta- 
tion or due to doping during epitaxy 0, 0, EH EH] ■ If the 
impurity concentration is high, a generation of nonequi- 
librium interstitial impurity atoms is possible within this 
layer during thermal treatment. These nonequilibrium 
interstitial atoms can diffuse before they transfer to the 
substitutional position or are trapped by immobile sinks. 
In many cases (see, for example, impurity profiles in 
@; EH EH EB|), the distribution of impurity atoms in 
an as-grown structure can be described by the Gaussian 
function 



C Q {x) = C(x, 0) = C m exp 



{x ~ Rpf 



2AiL : 



(27) 



where C m is the maximum value of impurity concentra- 
tion; R p is the position of the maximum, and Ai? p is the 



standard deviation. 

Let us suppose that the Gaussian function can also be 
used to describe the spatial distribution of the generation 
rate of impurity interstitials: 



G AI (x,t)=g m 1 exp 



AI . 



{x - R P f 



2AR> 



(28) 



where g m ! is the maximum rate of generation of intersti- 
tial impurity atoms. 

Taking into consideration expressions (|2"6")l and [)28p 
yields 



AI AI 

w{i,t) = % exp 

AI 



(A - Rp) 



2AiL 2 



S'(x B -0Cf. 



(29) 

Substituting the Green function (f25|) and the standard- 
izing function (|29|) into expression f9| allows one to ob- 
tain a spatial distribution of diffusing interstitial impu- 
rity atoms: 
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Calculating the integrals on the right-hand side of ex- 
pression (|30[) . one can obtain an explicit expression for 



the quasistationary distribution of interstitial impurity 
atoms: 



J 



exp(ui) 



d AI cosh vS + I, 
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where 
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(33) 
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(35) 
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Iai 



(42) 



(43) 



Postulating that the loss of substitutionally dissolved 
impurity atoms is equal in modulus to the rate of gener- 
ation of impurity interstitials (G AS (x.t) = — G AI (x, t)), 
taking into account that the distribution of impurity 
interstitial atoms C AI (x,t) = C AI (x) for the time- 
independent generation rate (|27p . and substituting ex- 
pressions (|2"7|) . (f2"8"|) , and (f3"Tj) into ©, one can calculate 
the concentration of impurity atoms in the substitutional 
position 



C(x,t) = —C AI (x) + C m (l- P /ll )exp 



AI\ 



(x - R P f 
2AR£ 



(44) 

where p AI is the fraction of the impurity atoms which 
transferred from the substitutional position into the in- 
terstitial one. 

Expressions (|3 1 1) and (|44]) are the obtained solution 
of the boundary value problem under consideration and 
can be used for verification of approximate numerical so- 
lutions and for simulation of interstitial diffusion. 
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(39) 



(40) 



(41) 



IV. SIMULATION 



The analytical solutions (|3 1 1) and (|44|) can be used for 
modeling different diffusion processes in semiconductor 
substrates. Below, we consider the case of nitrogen dif- 
fusion in gallium arsenide investigated in Nitrogen 
distributions of an as-grown and subsequently annealed 
specimen measured by secondary ion mass spectrometry 
(SIMS) are presented in Fig. [TJ In Q a GaAs layer with 
about 1 urn. thickness was grown on a gallium arsenide 
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substrate at 580 °C by solid source molecular beam epi- 
taxy. Intermediate introduction of an N 2 flow from a 
plasma source resulted in a buried N doping layer with 
a peak concentration of about 10 7 /tm -3 and a width of 
several 10 nm. This doping layer was thus sandwiched 
between a buffer and a cap layer of either roughly 0.5 /im 
thickness. Diffusion annealing was performed at a tem- 
perature of 822 °C for 15 h in sealed quartz ampoules. It 
can be seen from Fig. [1] that the nitrogen concentration 
profile after annealing is characterized by two extended 
low concentration "tails" directed into the bulk of the 
semiconductor and to its surface. 

Let us consider a possible mechanism of the forma- 
tion of such "tails". In recent years the mechanism of 
dopant diffusion in silicon crystals due to the forma- 
tion, migration, and dissociation of the "impurity atom 
- vacancy" or "impurity atom - self-interstitial" pairs 
(the pair diffusion mechanism) has become commonly 
accepted (see, for example [l|, [IE Qj], EE])- It is sup- 
posed within the framework of the pair diffusion mech- 
anism that a local thermodynamic equilibrium prevails 
between substitutionally dissolved dopant atoms, intrin- 
sic point defects, and the pairs. However, boron diffusion 
in silicon is often considered within the framework of the 
substitutional- interstitial mechanism [IE HE HH, when 
the silicon self-interstitial displaces an immobile impu- 
rity atom from the substitutional to the interstitial posi- 
tion. A migrating interstitial impurity atom in turn re- 
places the host atom becoming substitutional again (the 
so-called "kick-out mechanism"). 

It is supposed that the kick-out mechanism is also re- 
sponsible for the diffusion of gold in silic on M, m ; zinc 

m, HE HI HE HE HI HI, HI , nitrogen d Mam 

and beryllium 



27, \X 



jjnag- 
I, 1H in 



, HI, HI [2U2^j23, HI, Hi , nitrogenJU^ 
nesium |32l . 

GaAs. As follows from the investigations presented in 
[27l HE H?J , beryllium diffusion in other compound semi- 
conductors is governed by the kick-out mechanism too. 

To describe diffusion due to the kick-out mechanism it 
is also supposed that there is a local equilibrium between 
substitutionally dissolved dopant atoms, self-interstitials, 
and interstitial dopant atoms. In this case a mathemati- 
cal description of diffusion due to the kick-out mechanism 
is equal to the description of the pair diffusion [34[ . It was 
shown in [E HE H^l within the framework of diffusion in 
silicon governed by the "dopant atom - self-interstitial" 
pairs that the formation of the extended "tail" in the 
low concentration region of phosphorus profile occurs if 
the distribution of silicon self-interstitials is nonuniform, 
namely, the distribution of self-interstitials in the neutral 
charge state have to be nonuniform [5j, |38| . Thus, "tail" 
formation can be attributed to the nonuniform distribu- 
tion of self-interstitials in the neutral charge state. 

In Ref. it was supposed that nitrogen diffusion in 
GaAs is governed by the kick-out mechanism with inter- 
stitial As as a native point defect (mathematically equiv- 
alent to the diffusion due to the "impurity atom - As 
self-interstitial" pairs), and that two extended low con- 
centration "tails" on the nitrogen concentration profile 



are formed due to the nonuniform distribution of As in- 
terstitial atoms. In contrast to [3|, it is supposed in this 
paper that during thermal annealing the generation of 
nitrogen interstitials occurs in the buried layer and that 
migration of these nonequilibrium interstitial atoms is re- 
sponsible for the nitrogen redistribution. In our opinion, 
this mechanism of diffusion is more preferable because 
the atomic radius of nitrogen is smaller than arsenic ra- 
dius and hence in the doped GaAs the nonequilibrium 
nitrogen interstitials prevail rather than the arsenic ones. 
In addition, it was shown in [40| that the mass action law 
and local thermodynamic equilibrium between substitu- 
tionally dissolved dopant atoms, vacancies, and vacancy- 
impurity pairs are not valid in the low concentration re- 
gions of the abrupt dopant profile formed by low energy 
ion implantation. This conclusion is also true for dif- 
fusion due to the "dopant atom - self-interstitial" pairs 
and due to the kick-out mechanism. It is important to 
note that the buried nitrogen layer investigated in [3[ is 
very narrow and disequilibrium between substitutionally 
dissolved impurity atoms and diffusing species may well 
occur. Therefore, we try to explain the experimental data 
of [3| within the framework of migration of nonequilib- 
rium nitrogen interstitials. 

In Fig. Q] the results of simulation of nitrogen diffusion 
in GaAs obtained on the basis of analytical solution (f3Tj) 
and ([31]) are presented. 
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FIG. 1: Calculated nitrogen concentration profile (solid line) 
after thermal treatment of a GaAs substrate containing a 
buried nitrogen layer with a width of several 10 nm at a tem- 
perature of 822 °C for 15 h. The experimental data (dots) 
are taken from Stolwijk el at. [3]. 

The following values of simulation parameters were 
used to fit the calculated curve to the experimental nitro- 
gen profile: i) the parameters of the as-grown nitrogen 
distribution: C m = 0.66xl0 7 /mi~ 3 ; Rp = 0.468 /im; 
ARp = 0.016 /im; ii) the parameters of the nitrogen in- 
terstitial diffusion: l AI = 0.16 /mi; r A/ =0.1 s; g AI 
102.0 fim^s- 1 ; p AI = 0.84; vf ff = 1.5 /1ms" 1 . 

As can be seen from Fig.[TJ the calculated curve agrees 
well with the measured nitrogen concentration profile. 
Thus, the experimental data [3( can be explained on the 
basis of the migration of nonequilibrium nitrogen inter- 
stitials. It follows from the value of the fitting parameter 
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p AI that approximately 84% of the nitrogen atoms from 
the buried layer are being transferred to the transient 
interstitial positions. Migration of these noncquilibrium 
interstitial atoms results in the formation of two extended 
"tails" on the nitrogen concentration profile. 

V. CONCLUSIONS 

The analytical solution of the equations that describe 
impurity diffusion due to migration of noncquilibrium 
impurity intcrstitials is obtained for the case of Robin 
boundary condition on the surface of a semiconduc- 
tor. Using this solution, one can verify the correctness 
of the approximate numerical calculations obtained by 
the codes intended for simulation of diffusion processes. 
Moreover, it is possible to carry out an analytical simu- 
lation of a number of diffusion processes which are based 



on the migration of impurity interstitial atoms and can 
be used in the fabrication of semiconductor devices. The 
solution obtained can also be useful for simulation of im- 
purity diffusion in the doped layers with decananometer 
thickness because in these layers disequilibrium between 
immobile substitutional^/ dissolved impurity atoms, mi- 
grating self-interstitials, and migrating interstitial impu- 
rity atoms can take place. As an example, nitrogen redis- 
tribution from a buried layer during thermal annealing 
of GaAs substrate at a temperature of 822 °C for 15 h 
have been simulated. By comparison with the experi- 
mental data, the values of the parameters that describe 
interstitial nitrogen migration have been derived. For ex- 
ample, it was found that for the process under consider- 
ation approximately 84% of the nitrogen atoms occupied 
transient interstitial positions and the average migration 
length of these interstitial impurity atoms was 0.16 ^m. 
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